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1 Introduction 

It is well-known that a class of "almost integers" can be found using the theory of modular functions, 
and a few spectacular examples are given by Ramanujan [1]. They can be generated using some 
amazing properties of the j-function. Some of the numbers which are close approximations of 
integers are exp(7r-v/163) (sometimes known as Ramanujan's constant), exp(7r\/37) and exp(7r\/58)- 
These irrationals come close to an integer as follows: 



exp(7r\/37) = 199148648 - 0.219... • lO""^ 
exp(7r\/58) = 24591257752 - 0.177... • 10"^ 
exp(7r\/l63) = 262537412640768744 - 0.749... • 10^^^ 

Another surprising result comes from the average length of a segment in an isosceles right triangle 
with catheti of unit length. If / is this average length, then 

I = — (2 + 4\/2 + (4 + \/2) sinh^^(l)) = 0.4142933026... = {^2 - 1) - 0.8... • 10"^. 
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Such astonishing non-equalities are usually called almost identities or non-identities. Many exam- 
ples of such unexpected behaviour are known The four examples above are however different 
in essence: the first three come from a deep property of a complex mathematical object (the 
j-function) and the last has a good chance to be a genuine arithmetical coincidence. 

A natural question that comes to mind in presence of such a non-identity is therefore whether 
or not the phenomenon is purely coincidental, or comes from a more subtle process. For instance, 
in the equation 

e'^ -vr = 19.999099979... , 

it is not clear at all whether the almost identity pops up from a deep connection between e and tt 
or just because the expression happens to be close to 20. 

Recently, J.M. Borwein and P.B. Borwein discovered several families of almost identities 
leading to a systematic study of such phenomena. These were based on mathematical concepts 



that lead to clear explanations. Among the non-identities studied by these authors, let us mention 
the following striking example: 

oo ^ I — 

10(fc/ioo)2 - 100^ 1^^, 

correct to at least 18,000 digits. In this situation, the almost identity is not a coincidence. From 
the same viewpoint, let us mention as well the sequence 

n! 

" " 2(ln(2))"+i' 

for 1 ^ n ^ 17, discovered by D. Hickerson. These numbers are close to integers due to the fact that 
the above quotient is the dominant term in an infinite series whose sum is the number of possible 
outcomes of a race between n people (where ties are allowed). See [3] for the exact expression of 
these numbers. Here, once again, no coincidence. 



While we were studying the function 

oo 
k=l 

that appears in the analysis of the complexity of the binary gcd algorithm, we came to find a new 
family of almost identities. Let us define the real numbers n„ as follows: 

oo J 

The following equalities show the very strange behaviour of the almost identities generated by the 
sequence {un}- 



Ul = 


vr + 0.53... 


• 10-11 


U2 = 


1 + 0.48... 


• 10-10 


Us = 


5 + <^-^- 


. • 10-9 


U4 = 


\ + »«^- 


.10-9 


U5 = 


^ + 0.15.. 


. • 10-*^ 


Uq = 


1 + 0.29.. 


. • 10-*^ 



This article presents an explanation of this phenomenon and sheds light on the relation between 
Un and Un+2- We first study the cases with n = 1 and n = 2 by using the theory of Mellin transforms. 
From there, we exhibit the recurrence relation 

1 n-2\ 

7 7 Un-2 + Tn 

4 n — 1 J 
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and give the explicit values of r„ satisfying < ^ rio = 0.7227399... -10 , Vn G N. We also 
present a generalization of the phenomenon, leading to, e.g., the almost-identity 

oo -j^ 

M4)- E F^ = ^ + 0.82... -10-^ 

fc=— oo 

In this article, we will use the notation f{x) ~a g{x) for equivalent functions in a neighbourhood 
of a and log2 x for the logarithm in base 2 of x. Also, the set N is considered to contain the integer 
in the sequel. 



2 The cases n = 1 and n = 2 



The first two cases in our list are 



^i=ln(2). J2 



1 



A;=— oo 



2-k/2 + 2*^/2 



TT + 0.53... • 10-^^ 



and 



U2 = ln(2) • V rr 



1 + 0.48... • 10 



-10 



In the next section, wc will see that the expression of n > 2, can be explicited based on these 
first two almost identities. We therefore begin our study by these cases. Let us define the complex 
functions gi and 52 as 

gi{x) = —2 • (arctan(-\/x) — 7r/2) and g2(x) = — ^ — , > 

1 + X 

as well as the functions Gi and G2 defined as 



Gn{x) = J29ni'2''x) , sRx>0 , n = l,2. 



k=l 



The convergence of Gi is justified by the fact that in a neighbourhood of +00 we have 



poo pco / 

arctant — 7r/2 = — / ^^dv = — / ( 



1 



1 



2 A. ~^ a ~^ 
V V v" 



dv = 0{l/t). 



It ^ Jt 

The following equalities are justified because G\ and G2 converge uniformly on compact subsets of 
their domains, and therefore the derivative can be interchanged with the sum. Therefore, 



d 



hm [G„(2-«)] 



d 



lim T ^"(2^=-") 

a— »oo ' nil 



m.-^oo'- — ' du L' 
k=\ 



00 /'2(fe-m)y"/2 

lim ln(2) ■ y ^ ' 



lim ln(2) • y 



^ (^2-(fc-™)/2 + 2(fc-"^)/2)' 



(2.1) 
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where the hmit is understood with m E N. The game plan is then to express the functions Gi 
and G2 in a completely different manner in order to compute these limits. The keystone of this 
process is the Mellin transform Recall that the Mellin transform of a locally Lebesgue integrable 
function f{x) over ]0,cxd[ is the function 



ris) = / f{x)x'-ux. 

Jo 

The conditions f{x) ~o 0(x") and /(x) ~oo 0{x^), with u > v guarantee that f*{s) exists in 
the strip —u < ?R.s < —v. Mellin's inversion formula |H1 P-13] states that if / is continuous and 
c g] — u, —v[, then 

1 rc+ioo 

fix) = — / ns)x-'ds, 

and in a neighbourhood of 0, we have 

f{x)= J2 Res(r(s)x-^s^), 

where the summation is over the poles si of the function f*{s)x~^ whose real part is strictly smaller 
than c. 

Let g{x) be a locally Lebesgue integrable function over ]0, cxd[, /(x) = X^fcli and suppose 

that the convergence is uniform in ]0, oo[. Then 



/*(^) = / y29{2''x)x'-'dx 
k=i 

roo 



m. (2.2) 



Proposition 1 For x > 0, we have 



^ / \ ^ 1 / ^ / \ sin(2fc7rlog2(x)) 

= -2 - + - g fc. cosh (2^.7 ln(2)) 

where Si{x) is a power series in x, which converges in [0, 1[. 

Proof: As announced earlier, the idea is to use Mellin transforms in a back and forth process to 
reveal another expression of Gi . Using H2.2|) , we can write 

Gl{s) = (2.3) 
In order to compute gl, recall that in a neighbourhood of +00 we have arctani — 7r/2 = 0{l/t). 
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So, we can perform an integration by parts, as long as 3fts G ]0, l/2[: 

5*(s) = -2 / (arctan(^) - 7r/2) x^'^^da; 
Jo 

—2 • (arctan(-v/x) — 7r/2) • — 



-[^ /-oo ^s— 1/2 



2s 



l + x 



-dx 



S .Iq 1 + X 

IT 



-dx 



S cos ITS 

The last equality comes from the relation 



OO rf-S—l 



dx 



TT 



iQ 1 + X smvrs 
Using Mellin's inversion formula with c = 1/4 and H2.3() . we get 

-j^ ^l/4+ioo 



Gi{x) 



TT 



27rz ./i/4-ioo ^scosvrs/ 2^* — 1 



-ds 



Res 



TT 



X 



s COS vrs/ 2^* — 1 



-,si 



The poles of the function (j^jj^^) can be partitioned as follows: 

i) s = is a pole of order two, 

ii) the real simple poles —1/2 + /c. A; G Z, 

in) the imaginary simple poles Ikni/ ln{2), k £ Z \ {0}. 
The residues are then 

— 7rlog2(a;) — f at s = 0, 



(-2)''+2 



(l+2fc)(2'=+l-v/2) 

1 exp(— 2fc7ri log2(a;)) 
2l ■ fc-cosh(2fc7r^/ln(2)) 



V^a;'' at s = -1/2 - k , k e'M, 



at s = 2A:7rVln(2) , k G Z \ {0}, 



and the above sum becomes 
Gi{x) -- 



-| - 7rlog2(x) +^ — 
which proves the proposition. 



-2) 



k+2 



k=0 



(l + 2A:)(-2^'+i + ^2) 



'XX 



E 



sin(2A;7r log2(x)) 
^ A; • cosh(2fc7r2/ln(2)) 



□ 



Corollary 2 ui = vr + Xlfcli 



27r 



cosh(2fe7r^/ln(2)) ' 
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Proof: Based on (|'2.1() . we have 



Ul 



lim _[Gi(2-«) 



TT + lim 



27r 



cosh(2A;7r7ln(2)) 
and the last limit being equal to zero, the corollary is proven. 



□ 



The case n = 1 is then settled since the sum on the right-hand side of the equality of Corollary 
121 is in fact small: 



Ul 



2tt 



0.538914478... • 10 
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cosh (2A;7r2/ln(2)) 
Proposition 3 For x > 0, we have 

„ , s 1 , , , „ , , 27r ^ sin(2A;7r logoff)) 
G.{x) = --- log,(.) + S,ix) - ^ g ,inh (2fevrV?n(2)) 

where 5*2 (x) is a power series in x, converging in [0, 1[ such that 5*2 (x) = 0. 

Proof: The proof follows the same lines as in the first case. First, 

vr 



oo ^s-1 



-dx 



IQ \ + X SlUTTS 

and thus, once again based on (|2.2)) and (|2.3j) . we have 

»l/2+ioo 



G2{x) 



1/2-ioo 
1/2+ioo 

1/2-ioo 

Res 

5Rs,<l/2 



Gl{s)x-'ds 



IT \ X 



sinvrs/ 2'^ — 1 



-ds 



TT \ x 



V sinvrs/ 2^ — 1 

The poles of the function can be partitioned as follows: 
i) s = is a pole of order two, 
a) the real simple poles k, k ^'L \ {0}, 
Hi) the imaginary simple poles 2/i:vri/ ln(2), k G 'Z \ {0}. 
The residues are then 

— log2(x) — ^ at s = 0, 



_ {-2km 
1 ' In(2)-sinh(2fc7r2/ln(2)) 



^ . ^M-2kmlo^^{x)l_ ats = 2fcvri/ln(2), A:GZ\{0}. 
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The new expression of G2 is therefore 

k=l 



G,{x) = --- log,{x) - ^ - L ,inh (2A:vrV ln(2)) ' 



This concludes the proof. □ 

Corollary 4 u, = 1 + Zk=i sinh(2.lvin(2)) • 
Proof: We use here the same trick as in CoroUarylJJ 

U2 = lim -^[G2(2-") 



u=m 



— 1 r Iq { 

and the limit being equal to zero, the corollary is proven. □ 

Once again, this shows why the number U2 is almost an integer. Indeed the sum on the right- 
hand side is fairly small: 

27T °° 2kTr 
U2-l = r-—- V — , , , ^ . = 0.4885108992... • 10"^°. 
ln(2) ^ sinh(2/c7r2/ln(2)) 

3 The recurrence relation 

Having found the roots of the mystery related to the non-equalities ui 7^ vr and n2 7^ 1, we would 
now like to extend the method used in the previous section to understand why U3,U4, ... are so close 
to "good arithmetic numbers". Looking back to the cases n = 1,2, we see that the functions gi 
and (72 played a crucial role. The key was the fact that they satisfy the equalities 



d 



du 

The next lemma shows how we can extend them: 
Lemma 5 Let n e N, n > 2, and let 



In 2 

n = 1,2. 



Rn,k 



du, 



(2-(fe-«)/2 + 2(*^-")/2) 

1 / 2^^~^^ /"^ \ " ^ /l 



21n2- (n- 1) I 1 + 2^=-" / V 1 + 2'="" 

Then 

(I n-2\ 

in,k — 7 ■ 7 ln-2,k + Rn,k 

V 4 n — 1 ' 
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The proof is left to the reader, who can simply differentiate and check! The equality of the 
previous lemma can be used as follows. For n > 2, we have 



= ln(2)- lim V 



oo^ (2-{'=-H/2 + 2('=-H/2)" 



M2)- Hm2-[/„,.U 



m-+oo ■'— ' du 
k=l 



ln(2) • lim I y -r- 



\k=l 



1 n-2 



4 n-1 



In-2,k + Rn,k 



(3.1) 



Let us define, for n> 2, 

fn{x) -- 



'X \ 1 — X 



1+X 



1 + x 



, so that 



21n2- (n-1) 



/n(2 ") — -Rn,*;- 



If 



i^n(x)=^/n(2M, 



fe=l 



since this function converges uniformly and absolutely on compact subsets of ^x > 0, we can 
interchange derivation and summation to obtain 



d 



— -• lim -[F„(2-«)]^_^ = ln(2). lim {Y.:r^^''^l- 

2 • (n — 1) m-^oo du 'iu-m m-^oo \ f-^ du 



(3.2) 



Once again, wc use McUin transforms to find another expression for each of the functions Fn{x) 
in order to compute these limits. 

Proposition 6 The function F^, n > 2, can be represented as 



FJx) 



^ , , 47r sin(2/s7rlog2(a;)) , 

Snix) - V Cfc . , ^. , when n is even, 
ln2^ smh(2K7r"'/ln2) 

/-n, / N 47r , sin(2A;7rlog2(a;)) , . 
^Jxbn{x) — > Ofe — .1 ,m _9 It — when n is odd, 



In 2-^ ''cosh(2A;7r2/ln2) 



where Sn{x) is a power series converging in [0, 1[ such that 5(0) = 0. The coefficients Ck and hk 
are given by 



Ck 



n-=o(j^ + 4vr^feVln(2)^) 
(2/ -2)! 

27rfcn-=o((i + 1/2)^ + 47r^fcVln(2)- 
ln(2)(2Z - 1)! 



when n = 21, Z > 1, 



when n = 2Z + 1, I > 1. 
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Proof: First, 



as) 



fnix)x'' dx 

fix " 

+ 



CO n/2+s— 1 



in— 1 



dx. 



This expression can be evaluated with the help of the Gamma function T. Indeed, this function 
satisfies, see, e.g., PQ p. 47], 

mm ^ r 

r(p + q) Jo (1 + x)P+'i 

and therefore 

r(n/2 + s - l)r(n/2 - s) r(n/2 + s)T{n/2 - s - 1) 



-2s 



r(n-l) r(n-l) 
r(n/2 - 1 + s)r(n/2 - 1 - s) 

r(n- 1) ■ 



We used the equality r(n) = (n — l)r(n — 1) in the last step. Based on Euler's reflection formula 
r(s)r(l — s) = vr/ sin(7rs), see, e.g., p. 9], the previous equality leads to the following expressions, 
both correct for 5Rs g]0, 1/2[ : 



fm = { 



1-2 

vr 



Ylif - s^) when n = 2l, / > 1, 



{21 - 2)! sinvrs _^ 
i„ •^^n((j' + whenn = 2/ + l, Z ^ 1. 

(2/ — Ij! COSTTS 



The equality (|2.2j) and (|2.3|) lead once again to 

-1/4+joo 



/'i/4+JOO 

Fnix) = / F*(s)x-^ds 

Jl/4-ioo 

E Res(/;:(s)^,.,) 



5Rs;<l/4 

The poles of the function can be partitioned as follows: 

i) the imaginary simple poles 2/c7ri/ ln(2), k £ 'Z \ {0}, 

ii) the real simple poles: 

s = j €z Z with |j| > I — 2 when n = 21, 

s = j + 1/2 gZ+1/2 with \j\ > 1-2 when n = 2/ + 1. 
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Note that contrary to the cases we have seen so far, there are no poles with multiplicity. The 
real simple poles will clearly contribute to residues of the form a^x*^ when n is even and a^^/xx^ 
when n is odd. We do not exhibit the coefficients since we will not need them. The imaginary 
simple poles lead to residues at s = 2A;7ri/ln2 , A; 7^ 0, which are of the form 

^-2km/ In 2 



Res 



X 



2^* - 1 

The new expression of is therefore 

' ~^ .00 
47r 



,2A;7ri/ln2 = f*{2km/ \ri.2) 



Fn{x) 



E 

k=l-l 

00 



X ^ akx'' 



Esin(2A;7r log2(x)) 
... ^ sinh(2/c7r2/ln2) 

47r I, sin(2A;7r log2(x)) 



ln(2) ■ 

when n = 21, I > 1, 
when n = 2/ + 1, 1^1, 



ln2 ^ cosh(2A:7rVln2) 
k=l-l k=l ^ ' ' 

where the coefficients and 6fc are given in the proposition. This finishes the proof. 



□ 



Corollary 7 The sequence {u„}neN satisfies the following recurrence relation 

'1 n-2^ 



where 



2tt 



ln(2)(n 
2tt 



4 n - 1 

2kTT 



Un-2 + rn 



^'^^sinh(2A;7r2/ln2) 



ln(2)(n 



2k'K 



^ cosh(2A;7r2/ In 2) 



when n = 21, I > 1, 



when n = 2/ + 1, 1^1. 



Proof: Based on (|3.1j) . (|3.2|) and the previous proposition, we have 



1 n-2 



4 n - 1 



Un-2 



2 ■ in 



1 ) m-+oo du 



The limit in the above expression annihilates the limit of the power series of Fn and the only 
contributing term in the limit is the sinus series of This gives the expected expression of r„. 

□ 



The growth of the coefficients r„ is the combined effect of the increase of the values of and 
6/c and the decrease of (n — As a consequence, the sequence r„ is increasing for n ^ 10 and 

decreasing for n ^ 10, which gives 

< r„ rio = 0.7227399... • 10-^ 

We end this article by the following remark. The entire theory used here to explain why the 
numbers m„ are so close to elements in Q U vrQ has nothing to do with the presence of 2 in the 
denominator of 

1 
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One could argue that any sum of the type 

1 



ln(m) • 



fc=--oo 

has the potential to he close to Q or vrQ depending on the parity of n. As a matter of fact, we 
have, for example, 

oo ^ 

• E ^^^2^ = ^ + • lo-^ 

A;=— oo 

oo ^ 

1^(9) • E = vr + 0.15... • 10-^ 

A;=— oo 
oo ^ 

ln(4) • V ^ = 1 + 0.37... • 10"^ 

Based on what has been shown in this article, we can say that the "error" term is due to the size 
of ln(m) (in the hyperbolic functions of r„) and the smaller it is, the smaller the error will be. In 
other words, the choice m = 2 is the best one can do in order to maximize the resemblance with 
elements in O U vrO. 
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